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ABSTRACT 

We present a series of high resolution radio and optical observations of the CLASS 
gravitational lens system B1152+199 obtained with the Multi-Element Radio-Linked 
Interferometer Network (MERLIN), Very Long Baseline Array (VLBA) and Hubble 
Space Telescope (HST). Based on the milliarcsecond-scale substructure of the lensed 
radio components and precise optical astrometry for the lensing galaxy, we construct 
models for the system and place constraints on the galaxy mass profile. For a single 
galaxy model with surface mass density S(r) oc , we find that 0.95 < (3 < 1.21 at 
2a confidence. Including a second deflector to represent a possible satellite galaxy of 
the primary lens leads to slightly steeper mass profiles. 
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1 INTRODUCTION 

Multiple-image gravitational lens systems are potent tools 
for investigating a wide range of astrophysical and cosmolog- 
ical issues. One promising application of strong lensing is the 
determination of the Hubble parameter through the mea- 
surement of differential time delays (Refsdal 1964; Schechter 
et al. 1997; Kundic et al. 1997; Lovell et al. 1998; Wisotzki 
et al. 1998; Biggs et al. 1999; Fassnacht et al. 1999; Koop- 
mans et al. 2000; Patnaik & Narasimha 2001). In addition, 
lenses can directly constrain the inner several kiloparsecs of 
galaxy mass distributions (e.g., Kochanek 1991, 1995; Rusin 
& Ma 2001; Cohn et al. 2001; Murioz, Kochanek & Keeton 
2001). Much of the recent interest in finding arcsecond-scale 
gravitational lenses through systematic searches such as the 
Cosmic Lens All-Sky Survey (e.g., Myers et al. 1995, 1999) 
and the southern survey of Winn et al. (2000) has been 
driven by these goals. 

High resolution observations are the essential step in 
transforming lens systems from mere novelties into useful 
tools. For example, optical or near-infrared imaging with 
the Hubble Space Telescope (HST) is necessary to pinpoint 
the centre of the lensing galaxy, determine its shape and 
orientation, and search for nearby objects that may be per- 
turbing the potential. This is complemented by observations 
with high resolution radio arrays such as the Multi-Element 
Radio-Linked Interferometer Network (MERLIN) and the 
Very Long Baseline Array (VLBA), which offer precision 
astrometry and search for milliarcsecond-scale substructure 



in the lensed radio components. Taken together, these data 
provide the raw material for developing models of the grav- 
itational potential, from which the structure of the lensing 
galaxy is investigated and time delays are predicted for use 
in Hubble constant determination. 

B1152-I-199 (Myers et al. 1999) was identified as a grav- 
itational lens candidate in the third phase of the CLASS 
survey observations. The source consists of two compact ra- 
dio components with a flux density ratio of ~ 3:1, separated 
by l'.'56. Follow-up observations with the Very Large Ar- 
ray (VLA) demonstrated that the components have nearly 
identical two-point spectral indices between 8.4 and 15 GHz, 
as expected for images of a single lensed source. The lens- 
ing hypothesis was definitively confirmed via follow-up spec- 
troscopy obtained with the Keck II telescope, which detected 
lines from a background quasar at Zs — 1.019 and a fore- 
ground galaxy at Zd ~ 0.439. Imaging of the system with 
the Palomar 5-m telescope revealed a bright {nig — 16.5, 
rrii = 16.6 in Gunn magnitudes) stellar-like object at the ra- 
dio position. The lack of structure suggested that the fainter 
of the images may be attenuated by dust extinction in the 
lensing galaxy. Subsequent observations with the Nordic Op- 
tical Telescope (Toft, Hjorth & Burud 2000) detected a weak 
counterpart to the secondary lensed image. The large opti- 
cal flux ratio between the images (~ 60:1 at V band) points 
to signiflcant differential extinction by the lensing galaxy, 
which was detected following the subtraction of scaled point- 
spread functions (PSFs). 

This paper presents vastly improved radio and optical 
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observations of B1152+199, and constructs models for the 
lensing mass distribution based on the high resolution data. 
Section 2 analyses the radio properties of B1152+199, in- 
cluding expanded VLA radio spectra and deep observations 
with MERLIN and the VLBA. Section 3 presents optical 
imaging obtained with HST. In section 4 we investigate mod- 
els for the system and constrain the lensing galaxy mass 
profile. Section 5 summarises our findings and outlines ad- 
ditional work. 



2 RADIO IMAGING 

B 1152+ 199 consists of two radio components with a flux 
density ratio of ~ 3:1 and a separation of l'.'56 (Myers et al. 
1999). Each of the components is unresolved by the VLA. To 
investigate the spectral properties of the components over a 
wide range of frequencies, the system was observed at 1.4, 
5, 8.4 and 15 GHz using the VLA in A configuration on 
1999 July 15. The phase calibrator was J1150+242, and ob- 
servations of 3C286 were used to set the flux density scale. 
The VLA data sets were calibrated in the Astronomical Im- 
age Processing System (AIPS) using the standard procedure 
and analysed in DIFMAP (Shepherd 1997). In each case the 
visibility data were fitted to a pair of compact Gaussian com- 
ponents using several iterations of model-fitting and phase- 
only self-calibration, with a solution interval of 0.5 min. The 
component flux densities arc given in Table 1, and the ra- 
dio spectra are plotted in Figure 1. The spectra are nearly 
identical, as expected for lensed images of a single back- 
ground source. The overall spectral indices between 1.4 and 
15 GHz are al^4 = -0.07±0.01 (A) and af,^ = -0.07±0.01 
(B), respectively. The flatness of the radio spectra suggests 
that the source may exhibit sufficient variability to allow for 
the measurement of a differential time delay. However, trial 
monitoring observations using the VLA at 8.4 GHz have 
thus far failed to detect any significant variability (^ 2%) in 
the lensed components. 

MERLIN 5 GHz observations of B1152-I-199 were per- 
formed on 1999 January 2 for 7 hr and again on 1999 Jan- 
uary 5 for 14 hr. The flux density scale was determined by 

15 min observations of 3C286 and the point source calibra- 
tor OQ208. Alternate observations of the target source (8 
min) and a nearby phase reference source Jl 148-1- 186 (2 min) 
were iterated. The combined data were calibrated in AIPS 
and imaged in DIFMAP by repeating a cycle of cleaning 
and phasc-ordy self-calibration, starting with long solution 
intervals (40 min) and gradually decreasing to a minimum 
interval of 2 min. Once the model had sufficiently converged, 
an amplitude self-calibration was performed using a solu- 
tion interval of 30 min. The final map has an rms noise 
of 70 /iJy/beam and is shown in Figure 2. The data were 
modelled by two compact Gaussian components with flux 
densities of 52.6 ± 0.1 mJy (A) and 18.2 ± 0.1 mJy (B). No 
further emission was detected down to the 3(t level of the 
residual map. 

VLBA 5 GHz observations of B1152-I-199 were obtained 
on 2001 February 27 and March 18. The observing time 

was three hours per epoch, which was divided among three 
widely spaced hour angles to maximise the uv coverage. The 
data were recorded in two IFs, each of which was split into 

16 0.5-MHz channels giving a total bandwidth of 16 MHz. 
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Figure 1. B1152+199 component radio spectra based on the 
VLA 1.4, 5, 8.4 and 15 GHz data of 1999 July 15. 
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Figure 2. MERLIN 5 GHz observation of B1152+199 taken 1999 
January 2 and 5. The lowest contour is 3x the map rms noise of 
70 /LiJy/beam, and contour levels increase by factors of two. The 
restoring beam is 72.5 x 63.7 mas at +38.1°. The data have been 
naturally weighted. 



Table 1. B1152+199 component flux densities (in m,Iy) at 1.4, 
5, 8.4 and 15 GHz. Data for the 1999 July 15 VLA observation. 
Errors in the flux densities are taken to be equal to the rms noise 
of the respective maps. 



Comp 5i.4 
A 50.8 + 0.4 
B 17.1 + 0.4 



S5 SsA S\5 

55.2 + 0.1 48.5 + 0.1 43.0 + 0.3 
18.6 + 0.1 15.6 + 0.1 14.5 + 0.3 
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Figure 3. VLBA 5 GHz observation of B1152+199 taken 2001 February 27 and March 18. The lowest contour is 3x the map rms noise of 
75 /ijy/beam, and contour levels increase by factors of two. The beam is 3.6 X 1.9 mas at —1.8°. The data have been naturally weighted. 
Left: Component A. Right: Component B. 



Table 2. B1152+199 component substructure. Data from the 
VLBA 5 GHz observation of 2001 February 27 and 2001 March 
18. Listed are the coordinates of each model component, its flux 
density and model-fitted position angle (for the Gaussian jets). 
Errors on the positions are ~ 0.1 mas. All errors are the nominal 
1(7 uncertainties estimated by the AIPS OMFIT task. 



Comp 


Xl 


X2 


55 (mjy) 






= 


= 


33.4 ±0.2 




Ajot 


+0'.'0115 


-fo'.'oooe 


17.1 ±0.1 


+87° ± 1° 


Bcorc 


+0'.'9353 


-1'.'2454 


11.9 ±0.1 




Bjet 


+0'.'9394 


-1'.'2494 


5.0 ±0.1 


+143° ± 2° 



The data were digitized using two-bit sampling and both the 
Stokes L and R polarizations were recorded. The integration 
time per visibility was 2 seconds. A five minute scan of the 
bright source 4C39.25 was included for fringe-finding pur- 
poses. Alternate observations of the target source (3 min) 
and the nearby phase reference source J1148+186 (1.5 min) 
were then iterated. 

Data reduction was carried out in AIPS on the com- 
bined epochs. The data were first flagged and amplitude 
calibrated using recorded system temperatures and known 
telescope gain curves. Next the data were fringe-fitted, using 
4C39.25 to derive an initial delay correction. The phase cali- 
brator was then used to find the residual rate, phase and de- 
lay solutions as a function of time. This was initially fringe- 
fitted assuming that it was a point source and later mapped 
and self-calibrated to produce a model that included the 
effects of resolved structure. The solutions so found were 
then applied to B1152+199. The data were mapped using 
the IMAGR task. As the correlation centre of the data was 
a position close (50 mas) to component A, the data were 
not averaged in time or frequency so as to avoid any smear- 
ing in the final images. Maps were made of two fields at 
the positions of the lensed radio components. Several iter- 
ations of mapping and phase-only self-calibration were per- 
formed until the CLEAN model converged. The data were 
then amplitude self-calibrated using a long solution interval 



(80 min). The resulting maps are displayed in Figure 3 and 
have an rms noise of 75^Jy/beam. Each of the components 
is resolved into a compact core and extended jet structure. 
Finally, the data were model-fitted using the OMFIT task. 
For each image, two components were used to crudely de- 
scribe the observed structure: a delta function for the core 
and a Gaussian for the jet. The positions, fiux densities and 
orientations of these components are listed in Table 2. 

Note the apparent peculiarity of the B1152+192 ra- 
dio substructure. In component A the jet is quite straight. 
Specifically, the position angle of the Gaussian jet compo- 
nent (+87°) is nearly identical to that of an axis drawn 
between the centre of that component and the core. The jet 
in B, however, appears to be slightly bent. Here the position 
angle of the Gaussian jet (+143°) is significantly larger than 
that of an axis passing through the centres of the core and 
jet model components (+134°). The bending of B could be 
produced by substructure in the lensing mass distribution 
(e.g., Metcalf & Madau 2001). Deeper imaging at higher 
resolution will be essential for confirming this discrepancy 
and evaluating possible solutions. 



3 HUBBLE SPACE TELESCOPE IMAGING 

HST observations of B 1152+ 199 were obtained on 2000 
March 28 with the Wide Field and Planetary Camera 2 
(WFPC2). The F555W ("V-band") and F814W ("I-band") 
filters were used. For each filter, four dithered exposures of 
500 sec duration were obtained. A standard reduction was 
performed on the data sets using the Image Reduction and 
Analysis Facility (IRAF). The resulting I-band image is dis- 
played in Figure 4. The V-band image detected all the same 
features but at a lower signal-to-noise, and is not shown. The 
lensing galaxy (G) and counterparts to the lensed quasar 
components (A & B) are clearly detected. In addition, a 
faint extended emission feature (X) is seen to the west of 
G. Component X certainly cannot be an additional lensed 
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Figure 4. HST WFPC2 I-band image obtained 2000 March 28. 
The shading is logarithmic. Counterparts to the lensed quasar 
images (A & B), the lensing galaxy (G) and a possible satellite 
galaxy (X) are detected. The region displayed is 3" x 3". 



Table 3. HST photometry and astrometry. Listed are the posi- 
tions of the components derived from the I-band image, and their 
integrated magnitudes mv and mj . The magnitude of X was cal- 
culated using only the central 150 X 150 mas. Errors are 10 mas 
on the positions and ~ 0.3 on the magnitudes. 



Comp 


xi 


X2 


my 


mi 


A 


s 


= 


^ 17.3 


^16.7 


B 


-|-0^'942 


-1^'251 


21.9 


19.9 


G 


+0^'554 


-0'.'991 


22.6 


19.6 


X 


-0^'037 


-1'.'046 


25.0 


23.1 



image, as it has no radio counterpart. It may therefore be 
a satclhte galaxy of the primary lens. In addition, evidence 
of a diffuse arc is seen to the west of image B. This may be 
lensed extended emission from the quasar host galaxy. Its 
position would require that part of the host galaxy resides 
in the region enclosed by the tangential caustic curve, inside 
of which sources are lensed into four images. 

Photometry and relative astrometry were performed on 
the HST data. The positions and magnitudes of the com- 
ponents arc listed in Table 3. There were difficulties in 
subtracting PSFs because image A is probably saturated. 
Consequently, the errors in the magnitudes are quite large 
(~ 0.3). The relative optical and radio positions of A and B 
are consistent to within the HST error bars (~ 10 mas). The 
position of G differs from the ground-based measurement of 
Toft, Hjorth & Burud (2000) by more than 100 mas, which 
demonstrates the importance of HST observations for inves- 
tigating lens systems. Finally, elliptical isophotes were fit to 
G. At a distance of 225 mas from the galaxy centre we find 
a surface brightness axial ratio of / = 0.84 ± 0.9 and a po- 
sition angle of —63° ± 17°. Closer to the galaxy centre, the 



isophotes tend to have smaller axial ratios and somewhat 
more vertical orientations. 



4 MASS MODELLING 

We now investigate simple models for the B1152+199 lous- 
ing potential. The models were constrained using the posi- 
tions {xi, X2) of the two lensed core components from the 
VLBA data and the flux density ratio. The consistency of 
the measured flux density ratio among the various radio data 
sets strongly suggests that the value lias been affected very 
little by variability and is therefore close to the true magni- 
fication ratio. We thus set r = \Sa/Sb \ = 3.0 and assumed 
an uncertainty of Ar — 0.15 (5%) to account for possible 
effects from low-level variability or model-fitting errors. We 
set Axi = Ax2 = 0.1 mas, and the lens and source red- 
shifts to their measured values of 0.439 and 1.019, respec- 
tively (Myers et al. 1999). A flat Qa = 0.7 cosmology with 
Ho = 100/ikms-i Mpc-i was assumed for all calculations. 

Because the lensing mass in B1152-f 199 is dominated 
by a single galaxy, we modelled the system using a singular 
power- law ellipsoid (SPLE) mass distribution (e.g., Barkana 
1998), with scaled surface density 



k{xi,X2) = 



<1 



(xi + xiipyi^ 



(1) 



where /3 is the power-law slope, q is the normalisation and / 
is the projected axial ratio. Current constraints on the inner 
several kiloparcsecs of elliptical galaxies (e.g., Cohn et al. 
2001; Grogin & Narayan 1996; Kochanek 1995; Koopmans 
& Fassnacht 1999; Rix et al. 1997; Rusin & Ma 2001) favour 
mass distributions that are close to isothermal (/3 = 1), but 
a range of slopes (0.8 < /3 < 1.2) is also consistent with the 
data. We considered various values of /? in this regime. Our 
calculations made use of the rapidly converging scries so- 
lutions for the deflection angles and magnification matrices 
of power- law mass profiles derived by Barkana (1998) and 
implemented in the "FASTELL" software package. 

Lens modelling was performed with an image plane min- 
imisation (Kochanek 1991), which optimised the fit statistic 



X 



= E 



Ax? 



Ar2 



(2) 



where primed quantities are model-predicted and unprimed 
quantities arc observed. To reduce the rmmbcr of free pa- 
rameters and ensure a constrained model, the position of 
the galaxy was fixed at (+0'.'544, -0'.'991) relative to A. 
With the profile slope also fixed, the model has five pa- 
rameters: the normalisation, axial ratio and position angle 
of the SPLE, and two coordinates for the unlensed source. 
The number of degrees of freedom (NDF) is zero, so the 
best-fitting model should be described by = 0, regardless 
of the assumed uncertainties. Monte Carlo simulations were 
performed to examine the stability of the models. Gaussian- 
distributed errors accounting for observational uncertainties 
were added simultaneously to the galaxy position (10 mas 
in each coordinate) and flux density ratio (5%). The models 
were optimised for each of 10000 trials. The uncertainties in 
the parameters were found from the ranges enclosing 95% 
of the results. 

The optimised model parameters and time delays for 
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Figure 5. The critical curve (thick line) and caustics (thin lines) 
of the isothermal (/3 = 1) lens model. The filled circle marks the 
recovered source position. The open circles indicate the positions 
of the images. Dashed lines denote contours of constant time delay 
in increments of 9.0h~^ days outward from the global minimum 
at imago A. 

various values of /3 are listed in Table 4. The differential 
time delay for the isothermal (/3 = 1) model is 35.9 ± 2.0 

days in a fiat f^A = 0.7 cosmology. The delay varies 
linearly with the assumed profile slope, as predicted by Witt, 
Mao & Keeton (2001). Critical curves and caustics for the 
isothermal case are plotted in Figure 5. 

The above models can be subjected to a pair of con- 
sistency checks. First, the mass distributions of ellipticals 
tend to be well-aligned with the hght (Keeton, Kochanek & 
Falco 1998; Kochanek 2001), so it is useful to compare the 
position angles of the SPLEs with the surface brightness of 
the lensing galaxy. From Table 4 we see that the model- 
predicted position angles are generally compatible with the 
HST value, but the agreement diminishes as the profile is 
made steeper. The discrepancies are, however, broadly con- 
sistent with the range found by Kochanek (2001). Offsets 
between the mass and light position angles are often due to 
external shear fields, but we lack a sufficient number of con- 
straints to unravel the distribution of internal and external 
shear in B1152-I-199 at this time. 

Second, we considered the relative orientations of the 
B1152-I-199 radio components. An infinitesimal vector dxA 
associated with image A is mapped to a vector dxB associ- 
ated with image B according to 

dxB = Ag^ AAdxA (3) 

where Aa,b = dy/dx are the magnification matrices (e.g., 
Schneider, Ehlers & Falco 1992) at the image positions A 
and B, respectively, and y represents a vector in the source 
plane. We define the position angles of the radio components 
using the axes passing through the centres of their respective 
core and jet subcomponents, as our experience indicates that 



this will be the cleanest indicator of orientation. Given the 
values in Table 2, PAa = +87° ± 1° and PAb = 4-134° ± 
2°. Using eq. (4) and the recovered magnification matrices, 
we projected a vector dxA with PAa = -1-87° to a vector 
dxB and calculated its position angle. The resulting values 
are listed in Table 4, with the uncertainties determined by 
the Monte Carlo. Note that steeper mass profiles lead to 
larger position angles for component B. A mass profile with 
(3 ~ 1.1 provides the best fit to the observed orientations. 
The popular isothermal model is still compatible with the 
VLBA data, but shallower mass profiles fail to reproduce 
the substructure. 

To better quantify constraints on the mass profile, we 
calculated a new fit parameter 

2 2 . 2 . 2 .2 ( A\ 

X \img ,pos ' /Cimg ,flx ' /Cgal ,pos ' /CPAq 7 

that includes contributions from the position angle of com- 
ponent B (xpAb) and the galaxy centre {xlai,pos), which 
was now allowed to vary in order to best reproduce the sub- 
structure. The galaxy position was constrained by the HST 
astrometry, with a fit tolerance of 10 mas on each coordi- 
nate. The model-predicted position angle of B was again 
computed by fixing A at the observed value of -f87°. A fit 
tolerance of 2° was assumed for the orientation of B. The 
models have NDF = 1. A plot of Ax^ as a function of f3 
is presented in Figure 6. The value of Ax^ is dominated by 
the contribution from PAb ■ The data favour a lensing mass 
profile with 0.95 < /? < 1.21 at 2a confidence (Ax^ < 4). 
The limits are changed only slightly when the rather arbi- 
trary error assigned to the fiux density ratio is set to 1% or 
10%. Finally, the absence of additional components in the 
5 GHz MERLIN map requires that any third lensed image 
have a fiux density Sc such that Sa/Sc > 150 (see also 
Rusin & Ma 2001), where Sa is the fiux density of image A, 
and a conservative detection limit of 5 times the map rms 
noise level has been assumed. This condition is satisfied for 
all SPLE mass models with l3 > 0.60. Therefore, the lack 
of a detectable third image leads to no improvement in the 
profile constraints in this particular system, but is consistent 
with the bounds derived using the observed images. 

We also investigated the expected perturbation due to 
the possible satellite galaxy X. Relative velocity dispersions 
can be estimated from the galaxy luminosities L and Faber- 
Jackson (1976) relation a oc L^^*. Because q (x oc L^^^, 
qa /qx — 3 using the V-band data and go /qx — 5 using 
the I-band data.n We fixed the position of X according to 
the HST astrometry and modelled it as a spherical deflector 
with the same power-law slope as G. We ran trials with 
qc/qx = 3, 4 and 5, and the results are listed in Table 5. 
The inclusion of X has a negligible effect on the position 
angle and axial ratio of G for each /3, but the time delay is 
decreased by 15 — 25% relative to the single deflector case. In 
addition, the predicted value of PAb is smaller for a given 
/3, making steeper mass profiles necessary to reproduce the 
data in the two deflector model (Table 5). Fixing qc /qx ~ 4, 
for example, implies that 1.04 < 13 < 1.26 at 2cj confidence. 



* Note that these values assume spherical deflectors. Modifica- 
tions would be required in the presence of ellipticity (e.g., Keeton, 
Kochanek &; Seljak 1997). For this preliminary analysis, however, 
we ignore such issues. 
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Table 4. Lens modelling results for power-law deflectors with different mass profiles (3. Listed are the axis 
ratio /, galaxy position angle PAq (in degrees), predicted time delay At (in days), magnification matrix 
elements and Ag at the image positions, and the predicted position angle (in degrees) of component B 
assuming that PAa = +87° . Uncertainties are determined from the range enclosing 95% of the Monte Carlo 
results. All values assume a flat = 0.7 cosmology, and were calculated using the fit parameter defined in 
eq. (2). 

/3 = 0.8 ^ = 0.9 ^ = 1.0 13=1.1 P = 1.2 

f n7S7+°°2* 719+""'^"' 677+°°-*° 627+""'"' 

/ U-'*' -0.025 "^■'^''-0.030 ^•'J^'^-0.035 '-0.041 "•"^'-0.048 

PAg -TA.^i;l -^e.Sl^f, -79.4l|^ -82.6l4;4 -8.5.6l*;| 

At 28.3 ±f. 6 32.f±1.8 35.9 ± 2.0 ^^■'^^22 43.4+2 4 

A^,n +0.379l°°l? +0.4301°;°^^ +0.4821°;°^^ +0.5331^^, +0.5821°;°^^ 

Aa,12 -0.273too", -0.284±0.0f3 -0.290 ± 0.015 -0.292 ± 0.017 -0.292 ± 0.018 

Aa,22 +0.m2tilf^ +0.762l°°ii +0.838 ±0.013 +0.910l°-°J^ 

Ab,ii +0.247+°°^4 +0.317+°-°^° +0.3961°°^^ +0.484 + 0.073 +0.580l°°^^ 

Ab,,2 -0.68lt°;«f, -0.792l°;oi -0.905i°;°|? -1.0211^1? -1.137t°;°l« 

Ab,22 -0.356l»:<^;j -0.366t°;l?« -0.357l°;l^^r -0.324l°:l^^ '^^^ 

PAb +123.o1^;^ +126.i1^;3 +129.91*;^ +134.51*;^ +140.015;^ 




Profile slope /S 



Figure 6. Constraints on the lensing galaxy mass profile. Plotted 

is Ax'^ as a function of profile slope /3 for the single deflector 
model (solid line) and two deflector (dashed line) model with 
Ig/ix = 4. The horizontal marker indicates the 2cr confidence 
limit. The vertical marker denotes an additional upper bound for 
the two deflector model from the absence of a third image. 



For profiles steeper than isothermal, a third, negative-parity 
image can be formed near X. This image becomes brighter 
as /3 is increased over the range of profile slopes investigated. 
The absence of the image tightens the upper bound on the 
profile for more massive secondary deflectors: /3 < 1.13 for 
qa/lx = 3 and 3 < 1.24 for qc/qx = 4. 

Finally, note that if we had used a larger position an- 
gle for B, steeper mass profiles would be required to fit the 
data. Assuming that PAb = +143°, as derived from the 
orientation of the Gaussian model component in Tabic 2, 
the profile constraints become 1.14 < /3 < 1.34 (2a) for the 
single deflector case and 1.2 < /3 < 1.4 {2a) for the two de- 



Table 5. Lens modelling results for two deflector mass models. 
Listed are the fixed normalisation ratio ijg/<?X: the '2a confidence 
limit on the profile slope (/3 bound 1) from the observed images, 
an additional upper bound from the absence of a detectable third 
image (/3 bound II), and the time delay for the isothermal (/3 = 1) 
case in days. 



flector case with qc /qx =3 — 5. The absence of a detectable 
third image is satisfied for the same range of profile slopes 
given in Table 5. This rules out virtually all /3 for qc/qx =3 
and 4. 



5 DISCUSSION 

We have presented high resolution observations of the grav- 
itational lens CLASS B1152+199 obtained with MERLIN, 
the VLBA and HST. Such data are vital if lens systems are 
to fulfill their potential for constraining galaxy mass distri- 
butions or determining the Hubble parameter. While pre- 
vious spectroscopic (Myers et al. 1999) and ground-based 
optical (Toft, Hjorth &: Burud 2000) observations were suf- 
ficient to confirm B1152+199 as a gravitational lens, they 
offered little detailed knowledge about the structure of the 
system. The motivation for this paper was therefore to pro- 
vide the high resolution data necessary for more advanced 
astrophysical and cosmological investigations of this lens. 

VLA and MERLIN observations of B1152+199 show 
two compact radio components separated by l'.'56, with a 
flux density ratio of ~ 3:1. Both components have nearly 
identical flat radio spectra between 1.4 and 15 GHz, as 
expected for lenscd images of a single background source. 
VLBA 5 GHz observations resolved each of the components 
into a compact core and extended jet. HST observations 
with the WFPC2 camera detected optical counterparts to 
each of the quasar images. The flux density ratio is ~ 70:1 
at V band and ~ 20:1 at I band, confirming the strong ex- 
tinction of component B as it passes through the lens. In 



1g/<1x 13 bound I [3 bound II At 

3 1.06 < /3 < 1.26 13 < 1.13 27.7 

4 1.04 <I3< 1.26 /3 < 1.24 29.8 

5 1.03 < < 1.25 /3 < 1.32 30.8 
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addition to clearly detecting the primary lensing mass (G), 
the images show a faint nearby emission (X) feature that 
may be a satellite galaxy. 

We investigated mass models for B1152+199 based on 
the new radio and optical data. The predicted time delay for 
an ellipsoidal mass distribution with an isothermal profile is 
35.9 ± 2.0 days, assuming a flat SIa = 0.7 cosmology. 
The delay varies linearly with the assumed mass profile, as 
predicted by Witt, Mao & Keeton (2001). The inclusion of 
a second power-law deflector to represent the possible satel- 
lite galaxy X decreases the time delay by 15 — 25%. The 
relative orientations of the lensed radio components were 
shown to be powerful probes of the mass profile. Assuming 
a surface mass distribution E(r) oc r"'', the milliarcsecond- 
scale substructure requires that 0.95 < (3 < 1.21 {2a) for a 
single galaxy model, and slightly steeper profiles when X is 
included. The lower limit on /3 is particularly secure, as an 
alternate definition of the component orientations, based on 
the model-fitted position angles of the Gaussian jets, would 
lead to steeper mass profiles. B1152-fl99 is the simplest 
gravitational lens system in which both upper and lower 
bounds on the lensing mass profile have been derived. This 
system demonstrates the power of VLBI data for the inves- 
tigation of mass models. 

Additional modelling insights are likely to be obtained 
from higher resolution VLBI observations of B1152-I-199. 
Note that the jet associated with image A appears to have a 
pair of hotspots (Figure 3). If corresponding features can 
be resolved in component B, it may be possible to con- 
struct models based on the detailed substructure of the 
jets, similar to the analysis performed on QSO 0957-1-561 
(Barkana et al. 1999). Such observations should also deter- 
mine whether the slight bend in the B jet suggested by the 
current VLBA 5 GHz data is real. Observations with the 
HST Near Infrared Camera and Multi-Object Spectrometer 
(NICMOS) will also be undertaken in an attempt to detect 
and map emission from the quasar host galaxy. The informa- 
tion provided by lensed extended emission is a vital ingre- 
dient for properly reconstructing the mass distributions of 
lenses (Kochanek, Keeton & McLeod 2001). Finally, future 
studies will focus on measuring the differential time delay. 
The system is currently being monitored at the VLA along 
with several other new CLASS gravitational lens systems. 
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